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ABSTRACT 

The H atoms inside minihalos (i.e. halos with virial temperatures T v j r < 10 4 K, in the mass range 
roughly from 1O 4 M 0 to 1O 8 M 0 ) during the cosmic dark ages in a ACDM universe produce a redshifted 
background of collisionally-pumped 21-cm line radiation which can be seen in emission relative to the 
cosmic microwave background (CMB). Previously, we used semi-analytical calculations of the 21-cm 
signal from individual halos of different mass and redshift and the evolving mass function of minihalos 
to predict the mean brightness temperature of this 21-cm background and its angular fluctuations. 
Here we use high-resolution cosmological N-body and hydrodynamic simulations of structure formation 
at high redshift (z j> 8) to compute the mean brightness temperature of this background from both 
minihalos and the intergalactic medium (IGM) prior to the onset of Lya radiative pumping. We find 
that the 21-cm signal from gas in collapsed, virialized minihalos dominates over that from the diffuse 
shocked gas in the IGM. 

Subject headings: cosmology: theory — diffuse radiation — intergalactic medium large-scale struc¬ 
ture of universe — galaxies: formation - radio lines: galaxies 


1. INTRODUCTION 


One of the most promising means by which to observe 
the high redshift universe in the cosmic “dark ages” is 
through the 21-cm wavelength hyperfine transition of the 
neutral hydr ogen that is abun dant p rio r to rei oni zati on 
('e.g.lScott, Hr. Reei3l1990t iSuhra, mania,n Rr, Padmanahhanl 

1 19931) . Motivated by the prospect of new radio tele¬ 
scopes that will be able to observe such a signal, several 
specific observational techniques have been proposed. 
Amon g these are the an gular fluctu ations on th e sky 
e.g. IM adau. M eiksin. fc R.eesl [lfl^jTqzzj^f^ajHgOOnt 

”ii003; 


]ligY_£L_alj2i]fl2|- I SFM hereafter: iC iardi fc Mada 


Iliev et all 1 200. ‘{I iZaldarriaga. Furlanet to. fc Hernauistl 
2004 iRirlanetto. SokasiamfcHemauistl 12004) . features 
in the frequency spectrum of the s ignal average d over 
a substantial patch o f the sky l| Shaver et all 119991 
iGnedin fc Shaved 12004) and studies of absorption fea- 
tu res in the spectra of brigh t, high-re dshift radio sources 
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For most of these techniques, with the exception of 
foreground absorption against bright radio sources, the 
21-cm signal must be distinguished from the CMB, which 
is only possible if the 21-cm level population corresponds 
to a spin temperature Ts, which differs from the temper¬ 
ature, Tomb, of the CMB. Since radiative excitation and 
stimulated emission of this transition by CMB photons 
tends to drive the value of Ts toward Tcmb, some com¬ 
peting mechanism must exist to decouple Ts from Tomb- 
There are two main physical mechanisms by which the 
spin temperature is decoupled from the CMB tempera- 
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ture: “Lya pumping,” or absorption of radiation with a 
wavelength in the Lya transition, followed by decay into 
one of the hyperfine l evels of the ground state (the “F ield- 
Wouthuysen effect” - iWouthuvsenl 1 9521 IFieldl 1 9591) . and 
spin excha nge during collisions b etween neutral hydro¬ 
gen atoms (iPurcell fc F ield 1956). The efficiency of Lya 
pumping depends upon the intensity of the UV radiation 
field at the Lya transition, whereas the efficiency of col- 
lisional coupling depends upon the local gas density and 
temperature. 

For z ;> 150, these mechanisms are ineffective at de¬ 
coupling Tg from T 0 mb , since the kinetic temperature of 
the gas, Tk, is coupled to T 0 mb by inverse Compton scat¬ 
tering, and sources of UV radiation have not yet formed 
to initiate Lya pumping. At z 150, however, Tk drops 
below Tqmb and, for z ;> 20, gas at the mean density is 
sufficiently dense for collisions to couple Tg to Tk- Dur¬ 
ing the “dark ages,” therefore, when there is no Lya 
pumping, the mean 21-cm signal against the CMB will 
be zero at z ^ 150, then in absorption at 20 ^ z ^ 150. 

At lower redshift, collisions become negligible for gas at 
or below the cosmic mean density, and such gas becomes 
invisible until its spin temperature is again decoupled 
from the CMB by Lya pumping due to an early UV back¬ 
ground from the first stars and quasars. Even though 
collisional decoupling is ineffective for z <; 20 for gas at 
the mean density, gas in overdense and/or heated regions 
can still be collisionally-decoupled. In particular, the gas 
density within “minihalos” - virialized halos of dark and 
baryonic matter with masses 10 4 M 10 8 M 0 and 
virial temperatures T < 10 4 /i which are too low to col- 
lisionally ionize their H atoms - is sufficiently high so 
as to decouple its gas spin temperature from the CMB, 
with Tg > Tcmb in general, causing it to appear in 
emission (ISFM). ISFM predicted the mean and angu¬ 
lar fluctuations of the corresponding 21-cm signal by a 
semi-analytical calculation which integrated the equa¬ 
tion of transfer through individual minihalos of different 
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mass at different redshifts (z > 6) and summed these in¬ 
dividual halo contributions over the evolving statistical 
distr ibution of m inihalo masses in the ACDM universe. 
llliev et ]23 ij2nri3i extended these results to include non¬ 
linear biasing effects. These authors concluded that the 
fluctuations in intensity across the sky created by mini¬ 
halos were likely to be observable by the next genera¬ 
tion of low-frequency radio telescopes. Such observations 
could confirm the basic CDM paradigm and constrain 
the shape and amplitude of the power spectrum at much 
smaller scales than previously po s sible. 

Since then, iFurlanetto fc Loeb ( 120041 1 have suggested 
that shocked, overdense gas in the diffuse IGM (prior to 
the onset of Lya radiative pumping) is also capable of 
producing a 21-cm emission signal and that this IGM 
contribution to the mean signal will dominate over that 
from gas inside minihalos. Their conclusion is based 
on an extension of the Press-Schechter approximation 
(Press fc Schechter 1JZ4) that is used to determine the 
fraction of baryons in the diffuse IGM that are hot and 
dense enough to produce a 21-cm emission signal. We 
will address this question here. 

In order to quantify these effects, we have computed 
the 21-cm signal both from minihalos and the IGM at 
z )> 8 for the first time using high-resolution cosmolog¬ 
ical N-body and hydrodynamic simulations of structure 
formation. We have assumed a flat, ACDM cosmology 
with matter density parameter fl m — 0.27, cosmological 
constant Ha = 0.73, baryon density f h = 0.043, Hub¬ 
ble constant H = 70kms _1 Mpc _1 , linearly-extrapolated 
a 8h -i = 0.9 and the “untilted” Harrison-Zel’dovich pri¬ 
mordial power spectrum. 

In this paper, we present detailed, high-resolution gas 
and N-body simulations which predict the 21-cm signal 
at z > 6 due to collisional decoupling from the CMB 
before the UV background is strong enough to make de¬ 
coupling due to Lya pumping important. Because the 
Lya pumping efficiency is expected to fluctuate strongly 
until enough sources form to make the efficiency uniform 
(e.g. iBarkana fc Loebl 120041. the results presented here 
will also be relevant for isolated patches of the universe 
during reionization itself, which would depend upon the 
location and abundance of the first sources of UV ra¬ 
diation. Within such regions, we focus on properly re¬ 
solving the gasdynamics of structure formation at small 
scales through the use of high resolution gasdynamic and 
N-body simulations. We test the semi-analytical predic¬ 
tion of the halo model of ISFM for the contribution to the 
mean signal from gas in minihalos, and investigate the 
extent to which IGM gas may provide a non-negligible 
co ntribution to the total fluc tuating signal, as suggested 
bv IFurlanetto & Loebl (I2004H . 

These results were first summarized in lAhn et ID,-: 
(2!M). Here we shall describe our calculations in full 
and present our results in more detail. 

The outline of this paper is as follows. In §[2]we sum¬ 
marize the basic physics of the 21-cm emission and ab¬ 
sorption and the analytical model of ISFM. We also de¬ 
scribe our cosmological simulations and their initial con¬ 
ditions, and our method for obtaining the 21-cm signal 
from our simulations. In §|3we present our results. Our 
conclusions are summarized in §0] 

2. THE CALCULATION 


2.1. Physics of 21-cm signal from neutral hydrogen 

The hyperfine splitting of the ground state of hydro¬ 
gen leads to a transition with excitation temperature 
T* = 0.068 K, wavelength Ao = 21.16 cm, and frequency 
z/o = 1.417 GHz. The ratio of populations of the upper 
(ni) and lower (n o) states is characterized by the spin 
temperature Tg according to 

— = 3exp(-T*/T s ). (1) 


Neutral hydrogen at redshift z produces a differential 
signal relative to the CMB at redshifted wavelength 
21(l+z) cm only if Tg differs from Tcmb- The 21 cm tran¬ 
sition is seen in emission against the CMB if Tg > Tcmb 
or in absorption if Tg < Tcmb- The value of Tg is de¬ 
termined by the relative importance of collisional and 
radiative excitations. A hydrogen atom can 1) absorb a 
21-cm photon from the CMB (CMB pumping), 2) collide 
with another atom (collisional pumping) and 3) absorb 
a Lya photon to make a Lya transition, then decay to 
one of the hyperfine 21-cm levels (Lya pumping). These 
pumping mechanisms jointly determine the spin temper¬ 
ature, 


T S 


Tcmb + y a T a + y c T K 

1 + Va + Vc 


( 2 ) 


where T a is the color temperature of the Lyman-a pho¬ 
tons, Tk is the kinetic temperature, y a is the Lyman- 
a couplin g constant, and « r is the collisio nal coupling 
constant (IPurcell fc Fiddl 119561 |FiddH T9591. As seen in 
equation ©, the spin temperature deviates from Tcmb 
only when these couplings exist. Throughout this paper, 
we will consider only the collisionally coupled gas, or the 
case where y a = 0. This is valid when 1) light sources 
were not yet abundant enough to build up substantial 
Lyman-a background or 2) the region of interest is far 
enough away from light sources. The collisional coupling 
constant is given by 


_ x,c 10 

Vc T k A 10 ’ 


( 3 ) 


where Aiq = 2.85 x 10 _15 s _1 is the Einstein sponta¬ 
neous emission coefficient, and C io = k{1 — 0) nu is the 
atom -atom collisional de-excitation rateJPurcel^&Field| 
1956). We use k( 1 — 0) tabulated in lZvgelmanl (I2fl05f) 
which is valid for IK < Tk < 300K, and for higher T k 
we use k( 1 — 0) tabulated in lAllison & Dalgarno ^ Ii969fl 5 . 

The 21-cm line can be observed in either absorption or 
emission against the CMB, with a differential brightness 
temperature defined by 


5T,(zz) = Ti,(y) — Tcmb,o, (4) 

where r l\(v) is the brightness temperature at an observed 
frequency v and Xcmb.o is the present-day CMB temper¬ 
ature. T,(zz) satisfies the radiative transfer equation in 
the Rayleigh-Jeans limit, 


Tb(zz) = XcMB.oe 


dr„, 


' T s(^) 


1 + z' 


( 5 ) 


5 Definition of ni l— 0) in | All iscmfe Dal gar n Oj <T969h is not con¬ 
sistent with that in IZvgelmarJ ( 200^) • One should multiply their 
k( 1 — 0) b y (4/3) in order to calculat e C\q when using tabulated 
k(1 — 0) of Allison & Dalgarncl <196Sf) . 
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where r„ = fj" dr' v , is the 21 -cm optical depth of the 
neutral hydrogen atoms to photons in the CMB observed 

today at frequency 7 /, t' v , = J^"' dr",, is the optical depth 
of the neutral hydrogen atoms at redshift z' to photons 
at a frequency v' = v(\ + z') (the frequency which a 
comoving observer sees at redshift z'), Ts(z') is the spin 
temperature of intervening gas located at redshift z ', and 
the infinitesimal optical depth along the path of the pho¬ 
ton as it travels for cosmic time interval dt' is given by 


dr' v , =cdt' n v i(y',z') 
cdz' 


H{z')(l+z')\ |_ 32t tv' 


3c 2 Ai 0 uhi(2') 




T* 


Ta(z’)\ 


C6) 


where H(z) is the Hubble constant at redshift z, and nni 
is the local density of neutral hydrogen. The line profile 
4>(v r ) satisfies 


f* + oo 


dv'4>(y') = 1 , 


(7) 


and is in a general form. For instance, in the presence of 
both thermal Doppler broadening and the Doppler shift 
due to peculiar motion, the line profile is given by 







( 8 ) 


where A v' = (iz' 0 /c)^2kT/m, and the Doppler-shifted 
line center is given by 


v’ 0 = v o 


1 -\- (3 


(9) 


where /3 = v/c is the line-of-sight peculiar velocity of gas 
(in units of the speed of light) toward the observer. 


2 . 1 . 1 . The unperturbed universe 

Solutions to the general radiative transfer equation (eq. 
0) exist in simplified forms in limiting cases. If the line 
is un-broadened and un-shifted, i.e. </>(t/) = 5{y' — vq), 
the solution to equation © becomes 


Tb{v) = IcMB.oe 


Ts(z) 

1 + z 


1 — e 



( 10 ) 


where v and 2 satisfy vq = u(l + z), and 


r{z) s = 3 22 x 10 -= 


32nT s H(z) 




0.0224/ \0.135 


-0.5 


f Tqmb 


(i+^r 5 ,(n) 


using nni = xhi^h , where uh = 1.9 x 10~ 7 cm _3 (l + z) 3 
and xhi is the neutral fraction of hydrogen. The IGM ki¬ 
netic temperature, Tk, is coupled to Tqmb by Compton 
scattering at 2 ^ 134. For z <; 100, the kinetic temper¬ 
ature of the unperturbed gas in the ACDM universe is 
well approximated by 

T k ~T C mb{z = 134)(1 + z) 2 /(1 + 134) 2 

= 368.55K x (1 + z) 2 /{ 1 + 134) 2 . (12) 

Since Tcmb/ 7s < maxjl, T C mb/T k }, in general, t v < 1 
for the unperturbed IGM at all redshifts of interest (8 <; 



10 20 30 40 


z 

Fig. 1.— Analytical prediction for the mean 21-cm differen¬ 
tial brightness temperature due to collisionally-decoupled miniha¬ 
los and an unperturbed IGM. Shown are the results based on the 
Press-Schechter (solid) and the Sheth-Tormen (dotted) mass func¬ 
tions for halos and the contribution from IGM gas with cosmic 
mean density and temperature (dashed). In the bottom panel, 
we show the minihalo collapsed fraction in the ACDM universe, 
again based on the Press-Schechter (solid) and the Sheth-Tormen 
(dotted) mass functions. 


0 <; 50), as seen in equation (ITU) . In that case, equation 
m can be approximated as 


ST b (v) = [ T s {z ) - Tcmb(~)] t(z)/( 1 + z) (13) 


using equation 0 . We can use equations GH> and GZ) 
to describe the signal from the unperturbed gas in the 
Universe, because the line profile in this case is narrow 
enough to be approximated by a Dirac-<5 function. With 
equations ©, m and d, one can calculate 5Tb of 
the unperturbed gas of the Universe, well approximated 
by 


5T b {z) 


— 10 ml( 


TiogioU + z)] 

T[log 10 (l + 36.8)] ’ 


(14) 


where 


F(x) = dex(3172.36- 18037.3 x + 43430.3 x 2 - 57481.7 a ; 3 
+45150.0 x 4 - 21042.1 x 5 + 5388.50 x 6 - 585.165 x 7 ), (15) 


which is in absorption until collisional pum ping becomes 
neglig ible at z ~ 2 0 (Fig . UJ see also iBharadwai fc Alii 
1200 2 an cl IZvgelmanll200f?) . 


2 . 1 . 2 . Perturbed universe: optically thin case 

Thermal Doppler broadening and Doppler shift by pe¬ 
culiar motions would drive (f> to be broadened and shifted, 
causing overlap of line profiles. In such cases, the solu¬ 
tion to equation ©, in general, is not given in a simple 
form as in equation m■ We show here, however, that 
the simple solution given by equation m also applies 
to the overlapped line profile case, as long as the optical 
depths, both infinitesimal and integrated, of gas in the 
simulation box are small. In such optically-thin limit, 
equation © can be approximated as 

T b {v) = Tcmb,o(1 — T v) + f d T ' v , -7^—7- ( 16 ) 

J 0 1 + z 
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Fig. 2.— Map of the differential brightness temperature, 8Tb, (projected onto one surface of the box) for the redshifted 21-cm signal 
obtained from our highest resolution simulation, C4. Rows, top to bottom, show redshifts z=30, 20, and 10. Columns, left to right, 
represent contributions from minihalos, the IGM and the total signal. Note that the scale is linear in 8Tb f° r the upper two rows of images, 
but logarithmic for the bottom row. 


The differential brightness temperature from a simula¬ 
tion box at z with a redshift-spread A z (<C z) and angle- 
spread Afi is 


t \ — f dvdn6T b {is) 

5TlM = JTdP. 


(17) 


where the frequency and angle intervals of integration are 
set by the size of the box. As the angle integration is a 
simple sum of different line-of-sight contributions which 
do not interfere with each other, we can first perform the 
line-of-sight average, 




5T b {y) 


1.0.S. 


rlS-j- 

J V — 


v-\-8v/2 
8v/2 


dv 


(18) 


and then integrate over angles. In equation m , one can 
show that 


rv-\-8v j2 


dv = 


vo 


v-bv/l (1 + Z) 2 Jz-Az/2 


2 + A2/2 


dz\ 


(19) 


and 


rl'+ 8 v /2 pT u 

/ dv / c?r(,/ 

’ v — 8v j2 JO 


-2 


r A - /2 dZ ’^4^Fi, m 

z-Az/2 


H(z)(l + zr J z _ Az/2 32 ttT s (z') 


where we have used the fact that dv'(j)(v') = 1, and 
have also assumed that the thermal broadening and the 
Doppler shift by peculiar motions are negligible com¬ 
pared to the width of the box: (Ai/) t hermai •C Av^ ox , 
(Al/')peculiar < Azz( x)x . Using equations 1111, EH, and 
(EH), we then obtain 


6T b {v) 


1.0.S. 



T S(z') 
1 + 2 ' 


Tcmb.o ) t(z') 



( 21 ) 

which is simply an averaged superposition of contribu¬ 
tions of gas given by equation ra, along the line-of- 
sight. 

We have shown, in this section, that the 21-cm differ¬ 
ential brightness temperature ST b {v) can be calculated 
by a simple average of individual contributions from gas 
at different locations in a simulation box, as long as op¬ 
tical depths are small (eq. EH). Care needs to be taken, 
however, when gas achieves considerable optical depth. 
In § 12.1.31 we describe how one can calculate the sig- 


dz', 
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nal from optically-thick media, which are mostly located 
inside minihalos. 


2.1.3. Perturbed universe: minihalos 


Minihalos which start to emerge at z ~ 20 have tem¬ 
perature and density high enough to produce a significant 
emission signal (ISFM). As the optical depth through 
each minihalo is not negligible, the full radiative trans¬ 
fer equation (eq. 0) should be solved through individual 
minihalos. Once individual halo contribution Ai/ e fj5T b}l/0 
is obtained for each given halo mass M, one can calculate 
5Tb from all the minihalos at redshift z by integrating 
over the halo mass function dn/dM: 


_ r ( 1 max rj n 

5T b - C[ / Au e gST b „ 0 A —— dM , (22) 

voH(z) J Mmin dM 

where Av e s, 5T b ^ Vo , and A refer to the parameters 
of the individual virialized halo, Av a ft = Az/(1 + 
z ) _1 = (i'o/c) \/2kT/m{l + z) _1 , 5T bil/0 is the face- 
averaged differential brightness temperature at line cen¬ 
ter, and A is the projected surface area of the halo. 
ISFM based their calculation on the nonsingular, trun- 
cat ed iso thermal sph ere (TI S) m odel for CD M halos 
bv IShaniro. Iliev. fc Ragal (1199911 and llliev fc Shanird 
(2001), in which the halo density profile and virial tem¬ 
perature are fully specified by the halo mass and redshift, 
and the Press-Schechter (1974) mass function, which de¬ 
termines the number density of halos at a given redshift. 
The minimum minihalo mass M mnl is set by the Jeans 
mass 


Afj 


= 5.7xl0 3 


^ 0.15 


- 1/2 


n b h 2 \ 3/5 /i + x\ 

0.02 ) V 10 ) 


3/2 

Mq. 

(23) 


For M max , ISFM used the mass for which the virial tem¬ 
perature is 10 4 K: 


= 3.95X10' (2^) " ,2 (l±i) _3/2 M 0 . 

(24) 

The neutral baryonic fraction of halos with mass above 
M max is uncertain, because hydrogen will be partially 
ionized due to collisions and photoionization by internal 
sources. Thus, the mass range from M m ; n to M max natu¬ 
rally defines the mass range of minihalos which are fully 
neutral. Figure ^ depicts the predicted signals from un¬ 
perturbed gas as well as from minihalos. We show results 
for both t he Press-Schechter and the Sheth-Tormen mass 
functions dSheth fc Tormenll2002t) . 


2.2. Numerical Simulations 

We have run series of cosmological N-body and gas- 
dynamic simulations to derive the effect of gravitational 
collapse and the hydrodynamics on the predicted 21 cm 
signal from high redshift. Our computational box has a 
comoving size of O.S/i^Mpc, which is optimal for ade¬ 
quately resolving both the minihalos and the small-scale 
st ructure-form ation shocks. We used the code described 
in iRvu et all ill 99311 . which uses the particle-mesh (PM) 
scheme for calculating the gravity evolution and an Eu- 
lerian total variation diminishing (TVD) scheme for hy¬ 
drodynamics. We generated our initial conditions for 
the gas and dark matter distributions using the publicly 


available software COSMICS iMa fc Bertschingeill995i) . 
The N-body/hydro code uses an TV 3 grid and (TV/ 2) 3 
dark matter particles. In order to check the convergence 
of our results we ran simulations at different spatial reso¬ 
lutions, with grid sizes 128 3 , 256 3 , 512 3 , and 1024 3 , which 
we denote by Cl, C2, C3, and C4, respectively. We re¬ 
port our results in § 0 based on our highest-resolution 
sim ulati on C4 and discuss the convergence of the results 
in § 13.31 

After the decoupling of CMB photons from the bary¬ 
onic gas, the IGM gas cools adiabatically due to cosmic 
expansion (eq. El). Equation Cl a grees, fo r instance, 
with t he solution to the equation (1) in lBha.ra.dwa.i Hr. Alii 
(12001 which describes how T gas evolves exactly. This 
temperature, T gas (z), was used in the simulation to set 
the minimum temperature of baryonic gas, to avoid neg¬ 
ative temperatures 6 . If a gas cell is cooled below T sas (z), 
its temperature is set back to T gas (z). Such a temper¬ 
ature “floor” may overestimate the gas temperature of 
underdense regions, but because of their low density and 
temperature, y c is small in these regions. This implies 
that the spin temperature Tg would be very close to 
Tcmb; and their contribution to 5T b would also be neg¬ 
ligible, whether the kinetic temperature Xk is calculated 
accurately or not. 

In addition to the total 21-cm signal from our simu¬ 
lations, <52\igm, we are also interested in the relative 
contribution of the virialized minihalos and the IGM to 
the total signal, the sum of which gives the total 21-cm 
signal, 6T b , to t = 5T b , ha io + ST b , igm- First, we calculate 
the total mean signal as a simple average over the sim¬ 
ulation cells, 5T bit ot = J2i ^b,i/N 3 . The minihalo con¬ 
tribution is given by <5T bj ha io = fi^ T b,i/N 3 , where 
is the fraction of the DM mass in a cell i which is part of 
a halo. The IGM contribution can then be obtained as 

5T bt igm = 5T bt to t — ST bt haio = 5^(1 ~ fi)5T bti /N 3 . 

(25) 

In order to calculate the minihalo contribution to the 
total differential brightness temperature, JT^haio, one 
needs to first identify the halos in the simulation volume. 
We identified t he halos using a friends-of-friends (FOF) 
algorithm (iDavis et a,l.lll98ffl with a linking length pa¬ 
rameter of b = 0.25. The FOF algorithm applies to the 
dark matter N-body particles, rather than the gas in grid 
cells. Once this halo catalogue is processed for each time- 
slice of our N-body results, the baryonic component of 
each halo is identified for the grid cells of the hydrody¬ 
namics simulation which are contained within the volume 
of the halos in our FOF catalogue. We do this as follows. 
First, the density in each cell contributed by each DM 
particle is determined by the triangular-shaped cloud as¬ 
signment scheme. For each cell in which mass is con¬ 
tributed by the DM particles of a given halo, the gaseous 
baryonic component in that cell is assumed to contribute 
a fraction /, of its mass given by the fraction of the to¬ 
tal DM mass in that cell which is attributed to the halo 
DM particles. Accordingly, each cell i contributes an 
amount fi5T b ^ to the signal attributed to halo gas, while 

6 One should, in principle, use the locally varying minimum tem¬ 
perature. However, usage of a global minimum temperature is well 
justified as described in the text, and it is computationally cheaper 
than implementing a locally varying minimum temperature. 
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Fig. 3.— Volume-weighted probability distribution functions (PDFs) of gas density (l + <5) and differential brightness temperature ( 5T &) 
for C4 versus (l + (5) and 5T'b, respectively, for the total gas (solid), MH gas (dashed), and IGM gas (dotted). Also shown are the cumulative 
differential brightness temperatures, i.e. ( 5Tb(< [1 + <5])) and (5Tb(< 5Tb)). The top and bottom panels from left to right correspond to 
2=10, 20, and 30, respectively. 



Fig. 4.— Evolution of mean differential brightness temperature, 5Tb, of 21-cm background. (a)(left) Evolution of the total 21-cm signal 
vs. redshift. All data points are directly calculated from our highest resolution (C4) simulation box, with the assumption that optical depth 
is negligible throughout the box. (b)(right) 5Tb vs. redshift below 2 = 20. The contributions from minihalos (circles), the IGM (triangles), 
and the total (squares) are plotted, as labelled. For comparison, the result for the unperturbed IGM is also plotted (dashed-dot curves). 


(1 — fi)6Tb t i is assumed to be the signal from the IGM 
outside of the halo, where 5T^i is calculated from the cell 
as a whole. 

2.2.1. Semi-Analytical Calculation of the Halo 
Contribution 


Our numerical simulations have sufficiently high res¬ 
olution to find all halos in the computational box and 
the large-scale structure formation shocks, but not to re¬ 
solve the internal structure of the minihalos themselves. 
However, as ISFM have shown, in order to obtain the 
correct 21-cm signal from minihalos one needs to do a 
full radiative transfer calculation through each individ- 
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Fig. 5.— Semi-analytical minihalo signal vs. IGM signal. The 21-cm flux from each h alo in the simulation is found by modeling the 
internal structure and 21-cm line transfer through individual halos as described by ISFM (§ mu . to calculate the halo 21-cm signal from 
each halo more accurately. Same notation as in Figure^] The semi-analytical 21-cm minihalo emission is higher than the raw simulated 
minihalo signal in Figure[l] The IGM signal remains the same. The raw minihalo and total signals plotted in FigureQ](thin lines and open 
symbols) are also shown for comparison. 


ual minihalo density profile since, unlike the IGM gas, 
minihalos have a non-negligible optical depth at the 21- 
cm line. Hence, we can refine our estimate of the mini¬ 
halo contribution to the total 21-cm signal by combining 
our numerical halo catalogers with the semi-analytic cal¬ 
culation of individual minihalo contribution as found by 
ISFM. In their approach, as described in § ECU the gas 
den sity of each minihalo is assumed to follow a TIS pro¬ 
file (llliev & Shapirol200lD . radiative transfer calculations 
are performed to determine the 5Tb for different impact 
parameters, and, finally, the face-averaged 5Tb is calcu¬ 
lated (see ISFM, for details). The halo mass function, 
dn/dM , is provided by the halo catalogue we construct 
from the simulation. Each individual halo contribution, 
A v e ff5Tb t i, 0 A, depends on its mass and redshift of forma¬ 
tion (ISFM). Once we calculate Au e s6Tb tVo A, we then 
obtain the halo contribution using equation |22j . 

3. RESULTS 

3.1. Numerical 21-cm Brightness Temperature from 
Minihalos vs. IGM 

In this section, we describe the results from our simu¬ 
lations. In Figure0]we show (unfiltered) maps of the dif¬ 
ferential brightness temperature obtained directly from 
our numerical data for our highest-resolution simulation 
(C4), as described in § |3 We show the total signal, 
as well as the separate contributions from minihalos and 
IGM, derived as we described in 5 12.21 at redshifts z = 30, 
20, and 10. At z = 30, the earliest redshift shown (top 
row), most of the diffuse IGM gas is still in the quasi- 
linear regime and cold, thus largely in absorption against 
the CMB. At redshift z = 20 (middle row), the diffuse gas 
is still largely in absorption, while the (relatively few) ha¬ 
los that have already collapsed are strongly in emission. 
The combination of the two contributions creates a com¬ 
plex, patchy emission/absorption map, and absorption 
and emission partially cancel each other in the total mean 


signal. Finally, at z = 10 (bottom row), including its dif¬ 
fuse component, gas heated above Xcmb is widespread 
leading to a net emission against the CMB. The bulk of 
this 21-cm emission comes from the high-density knots 
and filaments. Although both the halo and IGM contri¬ 
butions come from roughly the same regions, the mini¬ 
halo emission is significantly more clustered, while the 
IGM emission is quite diffuse. 

In Figure^ we have plotted the volume-weighted prob¬ 
ability distribution functions (PDFs) for the gas density 
(1 + <5) and the differential brightness temperature con¬ 
tributions 5Tb as functions of each other. The PDFs for 
gas density show that, while the highest overdensities 
(5 )> 30) are typically found inside minihalos and the 
lower overdensities (5 <) 30) and underdensities (<5 <; 0) 
are typically associated with the IGM, there is some 
overlap of the distributions for these two components. 
A small fraction of the volume contains lower density 
minihalo gas and higher density IGM gas. However, the 
cumulative distributions show that these volumes hardly 
affect the total mean brightness temperature contributed 
by each component. Similarly, the PDFs for the bright¬ 
ness temperature show that, while the volume which con¬ 
tributes the highest brightness temperatures is predom¬ 
inantly inside minihalos and that which contributes the 
lower brightness temperatures is predominantly located 
in the IGM, there is, once again, some overlap of the 
PDFs. A small part of the IGM volume exhibits high 
brightness temperature, while a small part of the mini¬ 
halo volume shows low brightness temperature. Once 
again, however, the cumulative distributions show that 
these regions hardly affect the total mean brightness tem¬ 
peratures contributed by each component. 

In Figure 0] we quantify the relative contributions of 
the minihalos and diffuse IGM to the total mean 21-cm 
signal averaged over the whole computational box and 
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their evolution. The evolution roughly follows the naive 
analytical estimates, as was shown in Figure 0 The to¬ 
tal signal is deep in absorption, with 5Tb < —10 mK at 
z > 37. The 21-cm signal is completely dominated by 
the IGM contribution at this stage. The absorption sig¬ 
nal follows the analytical prediction for the unperturbed 
universe in S 12 . 1.11 well. since the density fluctuations are 
still small and the uniform-density assumption is reason¬ 
ably accurate. The absorption continually decreases as 
significant nonlinear structures start forming and por¬ 
tions of the gas became heated due to this structure for¬ 
mation. The net signal goes into emission after redshift 
z ~ 20, reaching up to ~ 5 mK by z ~ 8 . The emis¬ 
sion signal at 2 < 18 is due to both collapsed halos and 
the clumpy, hot IGM gas. In terms of their relative con¬ 
tributions, the minihalos dominate over the diffuse IGM 
at all times when the overall signal is in emission, be¬ 
low z = 18. We find _that the relative contributions 
to the total signal, \6T b j\ / (|£Ti,,haio| + |<U\igm|) 
where j means either “halo” or “IGM,” is nearly con¬ 
stant over two different redshift regimes: for z > 20 , 
|<57\igm| / ( 1 5Tb, halo | + | <U”mgm |) « 1 , while for z < 
16, |<57\ ha i 0 | / (|<57\haio| + |<5 T &jI gm|) ~ 0.7. In the 
transition region, 16 ^ z ^ 20 the relative contributions 
exhibit more complex behavior, approximately canceling 
each other out, resulting in a total signal which is close 
to zero. 

3.2. Refined Estimate of the Simulated Minihalo 21-cm 

Signal 

As we discussed in § 12 . 2.11 we can improve our purely 
numerical estimate of the minihalo 21 -cm signal by re¬ 
placing each halo’s flux with the value obtained by de¬ 
tailed radiative transfer calculations. We obtain the total 
minihalo signal from equation GU), with the theoretical 
mass function dn/dM replaced by the actual, numerical 
halo catalogue obtained from our simulations, and the 
individual minihalo contributions, Aiz e s5Tb^ 0 A, calcu¬ 
lated by modeling each halo as a TIS. 

We find that the resulting 21-cm signal from halos is 
stronger than the “raw” numerical signal obtained di¬ 
rectly from the simulated halos and dominates the over¬ 
all emission signal even more (Fig. 0 . This is despite 
the fact that our consideration of the more centrally- 
concentrated analytical density profiles increases the op¬ 
tical depth of each halo. We attribute this non-intuitive 
behavior to the fact that the density profiles of the mini¬ 
halos found in our simulations are not fully resolved. 
By modeling the halo density profiles in detail, the lo¬ 
cal density inside halos is boosted, which significantly 
increases the coupling constant y c , which, in turn, in¬ 
creases the total emission signal, even though the optical 
depth through each halo also increases simultaneously. 
Note that we use the same population of halos for both 
estimates, and only the internal halo properties are mod¬ 
ified. 

In FigureEl we show the total minihalo collapsed frac¬ 
tion obtained from the simulations compared to that 
from the theoretical PS and ST halo mass functions. We 
also show the minihalo contribution to the total differ¬ 
ential brightness temperature signal. We see that the 
collapsed fraction in minihalos in our simulation roughly 
agrees with the analytical predictions, mostly lying be- 
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Fig. 6.— A comparison of analytical and numerical minihalo re¬ 
sults. (a) (top) Differential brightness temperature of the 21-cm sig¬ 
nal from minihalos for semi-analytical model (dotted: with Press- 
Schechter mass function; dashed: with Sheth-Tormen mass func¬ 
tion) , simulat ion C4 numerical result (squares) and semi-analytical 
calculation f§ l2.2.H based upon simulation C4 mass function (tri¬ 
angles). (b) (bottom) Minihalo collapse fraction from simulation 
C4 (squares) and analytical mass functions (line types follow those 
of the top panel). 
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Fig. 7.— Numerical resolution convergence results. Mean differ¬ 
ential brightness temperature signals for simulations Cl (squares), 
C2 (triangles) and C3 (circles) in units of the corresponding signal 
obtained from our highest resolution simulation, C4. Shown are 
21-cm signals from the total (top panel), halos only (middle panel) 
and IGM only (bottom panel). 


tween the PS and ST results. On the other hand, the 
minihalo contribution to the total 21-cm background ob¬ 
tained directly from the simulation is below the theo¬ 
retical predictions based on either PS or ST mass func¬ 
tions. The agreement is restored, however, when we 
replace each minihalo contribution to the total flux by 
its analytically-modeled value. This, once again, un¬ 
derscores the importance of resolving the internal halo 
structure for correct predictions of their 21-cm emission. 
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3.3. Numerical Convergence 

We now compare cases Cl, C2, C3, and C4 to check the 
robustness of our results with respect to numerical con¬ 
vergence. In FigureO we show the differential brightness 
temperature signals for our three lower-resolution simu¬ 
lations, Cl, C2 and C3, in terms of the signal obtained 
from our highest resolution simulation, C4. We show the 
total signal, as well as each separate contribution, from 
either the halos or the IGM gas. At z > 20 most of the 
gas density fluctuations are still linear, and a change in 
the resolution barely affects the results. Thus a modest- 
resolution simulation, or even the analytical estimate for 
an unperturbed IGM, suffices to obtain reliable results. 
In contrast, at lower redshifts (z < 20) the results de¬ 
pend strongly on the resolution. The low-resolution sim¬ 
ulations Cl and C2 underestimate the resulting 21-cm 
signal significantly, by factors of up to a few. The re¬ 
sults from these simulations improve somewhat at lower 
redshifts, below 2 = 10, but results are still below the 
ones from C4 by ~ 30 — 50% and ~ 20% for simula¬ 
tions Cl and C2, respectively. This is true for either the 
minihalo, IGM or the total signal. The results from our 
medium-resolution simulation C3, on the other hand, are 
much closer to the ones from the high-resolution simula¬ 
tion C4, with the two generally agreeing to better than 
10%. This indicates convergence of our results to within 
a few per cent for simulation C4. Such behavior could be 
naively expected, since at z < 20 non-linear structures, 
both collapsed halos and mildly nonlinear, shocked IGM 
gas, form in abundance at the scales we are investigat¬ 
ing, and thus high resolution is required to resolve these 
properly, as our simulations confirm. 

The relative contributions of the minihalo and the IGM 
signals, on the other hand, show a more robust conver¬ 
gence. In all cases of different resolutions, we find that 
the minihalo signal dominates the IGM signal at z < 20, 
while the IGM signal dominates the minihalo signal at 
z > 20. For the purely numerical result, the relative 
contribution of minihalos to the emission signal is about 
70% at z < 15, peaks to 100% at z ~ 18, and drops to 
50% at z « 20. The exact value of the transition redshift 
varies slightly with resolution. For z < 14, minihalos 
contribute ~ 70% of the emission signal. For the case of 
semi-analytical calculation of the minihalo contribution 
based on the simulated halo catalogues, the relative con¬ 
tribution to the emission signal is slightly higher, ~ 75%. 

4. CONCLUSIONS 

We have run a set of cosmological N-body and hydro- 
dynamic simulations of the evolution of dark matter and 
baryonic gas at high redshift (6 < z < 100). With the 
assumption that radiative feedback effects from the first 
light sources are negligible, we calculated the mean dif¬ 
ferential brightness temperature of the redshifted 21-cm 
background at each redshift. The mean global signal is in 
absorption against the CMB above z ~ 20 and in overall 
emission below 2 ~ 18. At z > 20, the density fluctua¬ 
tions of the IGM gas are largely linear, and their absorp¬ 
tion signal is well approximated by the one that results 
from assuming uniform gas at the mean adiabatically- 
cooled IGM temperature. At z < 20, nonlinear struc¬ 
tures become common, both minihalos and clumpy, hot, 
mildly nonlinear IGM, resulting in an overall emission at 
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Fig. 8. — Numerical resolution convergence results. Relative 
contributions of minihalos and diffuse IGM gas to the total 21- 
cm background. The top panel shows the results obtained di¬ 
rectly from simulations (Cl: triangle, long-dashed; C2: square, 
short-dashed; C3: pentagon, dotted; C4: circle, solid). The bot- 
to m pan el shows the results which were semi-analytically refined 
(§ 12.2. II point- and line-types follow those of the top panel). 


21 -cm with differential brightness temperature of order 
a few mK. 

By identifying the halos in our simulations, we were 
able to separate and compare the relative contributions 
of the halos and the IGM gas to the total signal. We find 
that the emission from minihalos dominates over that 
from the IGM outside minihalos, for z 20. In par¬ 
ticular, the emission from minihalos contributes about 
70 — 75 % of the total emission signal at z < 17, peak¬ 
ing at 100% at z ~ 18, and balancing the absorption by 
the IGM gas at z ~ 20. In contrast, the absorption by 
cold IGM gas dominates the total signal for z > 20. 

These results appear to contradict the suggestion by 
IIArianetto"fc~LoeB (12004D . that the 21 -cm emission sig¬ 
nal would be dominated by the contribution of shocked 
gas in the diffuse IGM. They used the Press-Schechter 
formalism to estimate the fraction of the IGM outside of 
minihalos, which is shock-heated, by adopting a spher¬ 
ical infall model for the growth of density fluctuations 
and assuming that all gas inside the turn-around radius 
is shock-heated. This method is apparently not accurate 
enough to describe the filamentary nature of structure 
formation in the IGM. 

On the other hand, our results are consistent with the 
analytical estimates of the mean 21-cm emission signal 
from minihalos by ISFM. This indicates that the sta¬ 
tistical prediction of the collapsed and virialized regions 
identified as minihalos by the Press-Schechter formalism 
(or its refinement in terms of the ST formula), with virial 
temperatures T < 10 4 K, with halos characterized indi¬ 
vidually by the TIS model, is a reasonably good approx¬ 
imation for the mean 21-cm signal for minihalos at all 
redshifts and a good estimator even for the total mean 
signal including both minihalos and the diffuse IGM at 
z <) 20. This encourages us to believe that the angu¬ 
lar and spectral fluctuations in the 21-cm background 
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predicted by ISFM based on that model will also be 
borne out by future simulations involving a much larger 
volume than was simulated here. The current simula¬ 
tion volume is too small to be used to calculate the 
fluctuations in the 21-cm background because current 
plans for radio surveys to measure this background in¬ 
volve beams which will sample much larger angular scales 
(> arcminutes) than are subtended by our current box 
(A6»box ~ 0.2'(1 +^)§ 5 2 , where (1 + z) z > = (1 + z)/{1 + z')) 
and bandwidths (~ MHz) which are too large to resolve 
the depth of our simulation box in redshift-space (i.e. 

At/box ~ 40kHz(1 + z )ui 1 ^ 2 [T/(0.5fe~ 1 Mpc)1 ). Accord¬ 
ing to ISFM and llliev et alJ ( 2003 ). the fluctuations in 
the 21-cm background from minihalos are significantly 
enhanced by the fact that minihalos are biased relative 
to the total matter density fluctuations. A larger simu¬ 
lation volume than ours will also be necessary to sample 
this minihalo bias in a statistically meaningful way. This 
bias is likely to affect the minihalo contribution to the 
21 -cm background fluctuations substantially more than 
it does the diffuse IGM contribution, thereby boosting 
the relative importance of minihalos over the IGM even 
above the ratio of their contributions to the mean signal. 

We have considered the limit in which only collisional 
pumping is available to decouple the spin temperature 
from that of the CMB, and sources of radiative pumping 
have not yet emerged to compete with this process. 
The possibility exists, however, that an X-ray back¬ 
ground built up as sources formed inside some halos, 
which he ated the IG M while only partially ionizing 
it (e.g. lOh fc Haimanl 120031) . This heating might 
have boosted the kinetic temperature of the IGM and 
enh anced the e ffe ct o f collisional pumping there (e.g. 
iGhen fc Miralda-Escudel '200 li". Such X-ray heating 
would also have raised the minimum mass of minihalos 
which formed thereafter, filled with their fair share of 
neutral H atoms. When stellar sources began to form 
and build up the UV background at energies below the 
Lyman limit of hydrogen, Lya pumping could then have 
radiatively coupled Xjg to Ik, as well. The same sources 

7 Recently, IKnhlen. Mndan. K- Monty-ornery 1 2006 ) considered 
the X-rays emitted by an early miniquasar, finding that such an 
X-ray source can heat the IGM to as much as a few thousand 
degrees Kelvin without ionizing it. This boosts the 21-cm signal 
from collisionally-decoupled gas in the diffuse IGM significantly. 


presumably emitted UV radiation above the H Lyman 
limit, too, which ionized both the IGM and the minihalos 
within the HII regions s urrou nding these sources (e.g. 
[Sh^irmjlliey J _^jlaga|_j200^__|Il^w _Shai)iro. & R.agal 

l2005t llliev. Scannaoieco. fc Shaoirol I2005D . Such HII 

regions would have created holes in the 21-cm back¬ 
ground, which then originated only in the remaining 
neutral regions. Minihalos could have lost the neutral 
hydrogen gas responsible for their 21-cm emission, not 
only by “outside-in” photoionization by an external 
source, but also by “inside-out” ph otoionization by 
i ntern a l Pop III sta r formation (e . g. | Kit,a,va,ma, et all 
120041 lAlvarez. 'Brornm. Kr. Sha.nirnl 120061) . The H 2 
formation required for minihalos to form stars, however, 
is likely to have been suppressed easily by photodisso¬ 
ciation in the Lyman-Warner bands by the background 
of UV radiation created by the very first minihalos 
which formed stars, when the ionizing radiation back¬ 
ground was st ill much too l ow to cause reionization 
(iHaiman. Abel. V Reed , 2000 ). In that case, most 
minihalos would have remained intact until they were 
ionized from without. In the future, we plan to improve 
upon the current calculation by incorporating this more 
complicated physics. We also intend to run simulations 
with larger simulation boxes. This would allow us to 
predict the 21-cm fluctuation signal (e.g. ISFM) and 
determine whether the relative contribution of minihalos 
to the total signal, which we find to be about 70 75 % 

at z <; 20 for the mean signal, varies as the mean signal 
fluctuates. 
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